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Supersolidity in electron-hole bilayers with a large density imbal- 
ance 
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scopic, and nanoscale systems 
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density excitations) 
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Abstract. - We consider an electron-hole bilayer in the limit of extreme density imbalance, where 
a single particle in one layer interacts attractively with a Fermi liquid in the other parallel layer. 
Using an appropriate variational wave function for the dressed exciton, we provide strong evidence 
for the existence of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase in electron-hole bilayers 
with a large density imbalance. Furthermore, within this unusual limit of FFLO, we find that a 
dilute gas of minority particles forms excitons that condense into a two-dimensional "supersolid" . 



Introduction. — Pairing phenomena in two- 
component Fermi systems are a topic of fundamental 
interest, having relevance to a range of fields spanning 
superconductivity to QCD. For an attractive interspecies 
interaction, one can have Bardeen-Cooper-Schricffer 
(BCS) pairing in the weak-coupling limit or a Bose- 
Einstein condensate (BEC) of tightly-bound pairs in the 
strong-coupling limit. 

Of particular interest is the case where the densities 
of the two fermionic species are imbalanced, so that the 
interspecies pairing is then frustrated. Here, one ex- 
pects more exotic pairing scenarios such as the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) spatially-modulated 
phase [Ud], where fermions pair at finite centre-of-mass 
momentum, and both gauge invariance and translational 
invariance are spontaneously broken. However, an unam- 
biguous observation of the FFLO phase has remained elu- 
sive, despite being predicted more than four decades ago. 
The creation of spin-imbalanced Fermi gases in ultracold 
atomic gases OH] has recently revived the hope of real- 
ising the FFLO state, but, thus far, the atomic system 
has been dominated by phase separation between super- 
fluid and normal phases, with FFLO only occupying a 
tiny sliver of the predicted phase diagram in the three- 
dimensional (3D) case [5J[5]- 

On the other hand, electron-hole bilayers, where elec- 
trons and holes in a semiconductor are spatially separated 



into two closely-spaced quantum wells, may provide a bet- 
ter route to achieving the FFLO state. Here, electrons 
and holes can pair to form excitons which can then in 
principle condense 0[8]. Such bilayers have already been 
successfully produced by optically pumping coupled GaAs 
quantum wells (for recent experiments, see j9j[T0] ) . More 
recently, independently contacted layers have been fabri- 
cated in GaAs [HHH], where electron and hole layers can 
be separately loaded by biasing and doping, thus allow- 
ing the densities in each layer to be controlled individually 
and providing a means for generating a density imbalance. 
The additional advantage of considering these structures 
over optically-pumped coupled quantum wells lies in the 
extremely long lifetimes, where the exciton decay rate due 
to tunnelling recombination is essentially negligible and an 
equilibrium phase diagram can be regarded as accurate. 

The reason that electron-hole bilayers provide the ideal 
conditions for realising the FFLO state are twofold. 
Firstly, the reduced dimensionality of the bilayer system 
favours the FFLO phase over the normal phase owing to 
the enhanced Fermi-surface "nesting" [13], and secondly, 
the intra-layer Coulomb repulsion acts to suppress any 
macroscopic phase separation that might compete with 
the FFLO phase, in contrast to the cold-atom system 
where phase separation dominates. However, recent the- 
oretical work on density-imbalanced electron-hole bilay- 
ers |14H17j currently paints a rather uncertain picture of 
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FFLO. While some studies do predict regions of FFLO in 
the phase diagram, they either rely on an artificial tight- 
binding Hamiltonian with contact interaction |f 7] or on 
approximations that neglect screening |f 5j . Indeed, as we 
will show, screening crucially affects the stability of the 
FFLO phase. Other studies [HJQId] neglect the finite mo- 
mentum of the excitons in the FFLO phase and find that 
FFLO is often out-competed by other phases where the ex- 
cess particles simply coexist with the excitonic supcrfluid. 
Moreover, none of them consider competing phases such 
as charge-density wave (CDW) or Wigner crystal phases, 
which have been shown to be dramatically enhanced in 
mass-asymmetric electron- hole bilayers [18]. Clearly, fur- 
ther work is needed to ascertain the existence of FFLO. 

In this letter, we provide strong evidence for the exis- 
tence of the FFLO phase in an electron-hole bilayer with 
a large density imbalance. To address the problem in a 
more controlled manner, we consider the limit of extreme 
imbalance where we essentially have a single particle in 
one layer interacting attractively with a Fermi liquid in 
the other parallel layer. This allows us to rule out CDW or 
Wigner crystal phases induced by the presence of the other 
layer. It also enables us to include interactions between 
excitons that go beyond mean-field theory. From this anal- 
ysis, we expose an unusual bosonic limit of FFLO, where a 
dilute gas of excitons forms a condensate with a 2D spatial 
modulation, a phase otherwise known as a supersolid. Our 
strong-coupling version of FFLO should also be of interest 
to the communities working on strongly-correlated super- 
conductors, where one often encounters spatial textures 
coexisting with superconductivity. 

Model. — The basic Hamiltonian for a 2D electron- 
hole bilayer, H — Ho + Hi nt , consists of a kinetic part, 
Hq = X^kcr £ ko-cj co .c k(T , and an interaction part (fi is the 
system area): 
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Here, the labels a = {1,2} denote the different 2D lay- 
ers, and wc approximate the dispersions as quadratic, 
e ker = ft 2 k 2 /2m CT , which is reasonable for sufficiently small 
momenta in GaAs quantum wells. The bare Coulomb 
inter- and intra-layer interactions are respectively given 

by 
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where d is the bilayer distance, e is the electron charge, 
and e is the material dielectric constant. By introduc- 
ing twice the reduced mass m = 2(l/mi + l/m 2 )~ 1 , we 
can define the exciton Bohr radius ao = eh 2 /me 2 and the 
exciton Rydberg Eq = e 2 /ea . We also ignore the spin 



degrees of freedom and assume the layers are completely 
spin-polarised by a parallel magnetic field, but we shall 
return to this point later. 

We focus on the problem of a single particle in the 2nd 
(er = 2) layer interacting attractively with a Fermi liquid 
in the 1st (a = I) layer. This minority particle can either 
be a hole interacting with the electron layer or an electron 
interacting with the hole layer — the two cases are simply 
obtained by inverting the mass ratio a = mi/m\. For 
a fixed bilayer distance d/ao, the only other relevant pa- 
rameter is the dimensionless density r s = 2/fc^ao, where 
kp = 2y / 7m7 is the Fermi wave vector of the filled 1st 
layer with density n\. Note that this work focuses on es- 
tablishing the equilibrium phase diagram of a fully unbal- 
anced electron-hole bilayer, a problem somewhat different 
from the phenomenology of the x-ray-edge singularity |19j , 
where a dynamical transition is caused by a sudden local 
perturbation, e.g. the emission or absorption of a photon. 
Such a dynamical transition would also require a term in 
the Hamiltonian that transfers electrons between layers, 
which is clearly absent in (JXJ) . 

In order to determine the phase diagram of the fully 
unbalanced electron-hole bilayer, we need to establish the 
system ground state. To this end, we consider the fol- 
lowing variational state for an excitonic quasi-particle (or 
electron- hole pair): 
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where \FS) represents the Fermi sea of 1-particlcs filled 
up to wave vector hp, and we use the notation X)k>fe F = 
^fe>fe F ^ According to its definition, the excitonic wave 
function (pu.Q has relative momentum k and centrc-of-mass 
momentum Q. The spread in relative momentum k of ^kQ 
is set by the inverse of the size of the exciton bound state. 
Note that Eq. (J3j) gives a good description of both low and 
high density limits: In the low density limit, the state ([3]) 
coincides with the exact two-body (molecular excitonic) 
state in the vacuum limit (r s — > oo), while in the high 
density limit, where the particles barely interact, it gives 
the non-interacting "unbound" state: 
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corresponding to </?kQ = ^kfe F k^Qfe F k- ^ n tn is case, 
the minority particle is not bound to a majority particle 
and instead occupies a scattering state with well-defined 
(zero) momentum. We handle the regime between these 
two limits by considering screened interactions within the 
Random Phase Approximation (RPA). For the single- 
minority-particle case, the screened interactions U^ c , 
have the simple expressions 
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where the static polarisation operator ITi for the 1st layer 
is given by the Lindhard function: 
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with the number of particle flavours N s = 1 for the spin- 
polarised case. Considering screened interactions within 
RPA corresponds to effectively "dressing" the particles in 
our wave function with density fluctuations (i.e., an infi- 
nite number of particle- hole pairs). In addition, RPA pro- 
vides a good approximation in the long- wavelength limit, 
q <C 2kp, and therefore gives a good estimate of the un- 
binding transition: Here, as we will show later, the wave 
function ip^Q in (J3]) is strongly peaked at k = kp and so 
the momentum transfer is small. More generally, RPA 
should be reliable for sufficiently large d for the pairing 
coupling <7q C , since this removes the short-range interac- 
tions between electrons and holes where RPA has prob- 
lems [20]. 

The exciton ^(Q)) and unbound states resem- 

ble, respectively, the molecule and polaron wave functions 
used in ultracold atomic Fermi gases |21) . But there is an 
important difference: In cold atoms, one can accurately 
model the phase diagram using a wave function ansatz 
that explicitly includes just one particle-hole excitation 
on top of the non-interacting approximation to the Fermi 
sea \FS) for both the molecule and polaron wave func- 
tions [21][22]. This simple treatment is possible for cold 
atoms because the majority Fermi sea is non-interacting 
and the interspecies interaction is short-range. However, 
for long-range interactions and an interacting majority 
Fermi sea, this perturbative expansion is bound to fail, 
as it is never profitable to excite just one particle-hole 
pair. Indeed, the singular nature of the Coulomb interac- 
tion will generate an infinite number of particle-hole ex- 
citations at the Fermi surface. Thus, instead of explicitly 
including a small number of particle-hole excitations in 
the wave functions, we implicitly include an infinite num- 
ber of particle-hole excitations by replacing the Coulomb 
potentials in Eq. ([7]) with screened potentials and using, 
implicitly, an interacting Fermi sea. 

Phase diagram. — By minimising the expectation 
value (*(Q)| {H - E) |*(Q)) with respect to the ampli- 
tude iy9kq , we obtain an eigenvalue equation for the exciton 
energy E: 



5 £ V "A 

k'<fc F / 



\ J2 flklk'Pk'Q ( 7 ) 



k'>fc F 



Similarly to what was done in ultracold polarised Fermi 
gases [2TJ[22], we compare the energy of the excitonic 
molecular wave function |^(Q)) with that of the un- 
bound state )i i- e -i the exciton is bound if the low- 



est eigenvalue of ([7J is such that E < h 2 k F /2mi — 

H Sk'<fc F ^fe F k-k'- 

It is useful to note that the eigenvalue equation ([7]) co- 
incides, in the limit of full imbalance, with the mean- 
field gap equation employed to describe the BEC-BCS 
crossover in imbalanced electron- hole bilayers (see, e.g., 
Rcfs. P3HS])- To see this, we neglect the intra-laycr 
Coulomb repulsion and rewrite Eq. ([Jj in terms of the 
"gap", A kQ = i Ek'>fe F 5k-k'Vk'Q, therefore obtaining 
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This corresponds to the linearised version of the mean- 
field gap equation of Refs. [TJHS], which is just the form 
one would expect in this limit since the gap A^q be- 
comes macroscropically small as we approach full imbal- 
ance. To complete the correspondence between Eq. © 
and the linearised mean-field gap equation, we require 
that the chemical potential of the minority particles be 
\x% = E — h 2 k F /2mi, since the chemical potential of the 
majority particles is fixed to \i\ = Tck F j2m\. Thus, the 
conditions for a bound state outlined above require that 
fi2 < 0, implying that /X2 gives the exciton binding en- 
ergy in this limit, which also matches with the mean- field 
theory. 

Here, as for the excitonic wave function </?kQ, the or- 
der parameter A k Q describes a pair with centre-of-mass 
momentum Q. Therefore, in this context, it is natural to 
identify an exciton with minimum energy at non-zero mo- 
mentum Q with the FFLO phase in the large-imbalance 
limit — this will be justified further later when we con- 
sider the case of a dilute gas of minority particles. We 
emphasise that in this work, similarly to Ref. [15], and 
contrary to Refs. [T4lll6| . where only the limit Q — > has 
been considered, we allow the pair centre-of-mass momen- 
tum to be finite and minimise the energy with respect to 

Q 

The phase diagram for the fully-imbalanced bilayer sys- 
tem is obtained by converting the eigenvalue equation (J7J| 
into a matrix equation and numerically solving for the low- 
est energy eigenvalue. We fix the electron-hole mass ratio 
to 4, which is approximately its value in GaAs experi- 
ments, so that the relevant mass ratios are a = {0.25,4}. 
In the low-density limit, r s — > oo, we recover the two- 
body limit, and thus we expect a bound exciton with 
Q = (which we refer to as SF) for sufficiently large 
r s , as shown in Fig. Q] In the opposite limit, where r s 
is small, we see that the the screened interactions cause 
the exciton to eventually unbind and enter the "normal" 
(N) phase, as expected. However, the crucial point is that 
in a significant region of the phase diagram at interme- 
diate densities, the system ground state is a bound exci- 
ton with finite momentum Q, which we label as FFLO, 
as previously explained. The size of the FFLO region is 
greatest if the minority particle is an electron (a = 0.25) 
rather than a hole (a = 4). Indeed, the FFLO region is 
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Fig. 2: Behaviour of the momentum Qmin that minimises the 
exciton's energy as a function of r s for fixed d and a — m^/mi. 
The dashed line corresponds to Q = fc^ = 2/r s ao- Within the 
FFLO phase, Q m i n evolves continuously from Q = to Q = fop 
with decreasing r s . In the normal phase, Qmin = &f because 
the unbound state corresponds to tfikQ = <5 k fcj _, k <5q ;_ Fk . 



Fig. 1: Single- minority-particle phase diagrams as a function 
of the interaction parameter r s . In the top two panels, the 
bilayer distance d/ao varies and the mass ratio a = m^jmx is 
fixed to typical values in GaAs (a = 0.25, 4), while the bottom 
panel has a varying and d/ao — 1. In all cases, the inter- and 
infra-layer interactions have been screened using RPA. The 
superfluid (SF) region corresponds to excitons with centre-of- 
mass momentum Q = 0, while the "FFLO" excitons have their 
lowest energy when Q / 0. The normal (N) region is where 
there are no bound excitons. Refer to Fig [2] for the behaviour 
of Q within the FFLO region for d/ao — 1- 

generally enhanced (and shifted to larger r s ) when the mi- 
nority particle is lighter, as demonstrated in the bottom 
panel of Fig. [TJ The reason for this is simple: an exciton 
with Q = requires the minority particle to sit above the 
Fermi sea, but a small a increases the kinetic energy cost 
for this, thus favouring the formation of an FFLO exci- 
ton, where the minority particle can sit below the Fermi 
surface. 

The FFLO region is also enlarged by increasing the bi- 
layer distance d (Fig. [IJ. For large d, large momentum 
scattering |k — k'| > 1/d is suppressed in <7£1 k / and this 
tends to favour FFLO, where the wave function ipu.Q is 
peaked in the direction of Q, over SF. However, larger 
d also demands larger r s to achieve FFLO, and so we 
eventually expect Wigner crystallisation in the 1st layer 
to destroy FFLO. According to estimates from Quantum 
Monte Carlo (QMC) calculations [251121] . Wigner crystal- 
lization occurs once r s > 70a/ (1 + a). Thus, the distance 
d required to see FFLO sensitively depends on a, e.g. for 



a = 0.25, we ideally want d/ao < 1. 

Finally, we note that both the SF-FFLO and FFLO-N 
transitions are second order, with the exciton momentum 
Q varying continuously (Fig. [2]). Also, we always find that 
Q = kp at the FFLO-N transition. 

Interestingly enough, we find that the size of the FFLO 
region actually strongly depends on the screening of the 
interaction. In order to understand this, we consider the 
gap-equation form of the eigenvalue equation ^ for the 
case of unscreened interactions. After converting sums 
into integrals, and rescaling momenta by Oq 1 and energies 
by Eq, we obtain: 

= f d*V e-^'\ A k , Q 
kQ A<>£ 2tt |k-k'| -£ + e Q _ k ,, 2 + e k)>1 () 

Now, at the unbinding transition [E = Aa/r 2 s (l + a) in 
rescaled units] , we see that that the integral is logarithmi- 
cally divergent for Q = k' = k = 2k/r s and so we must 
take r s = for the equation to be satisfied. This implies 
that, for the bare Coulomb inter-layer interaction g q , the 
exciton with momentum Q = 2/r s = kp for r s C 1 is 
always bound, a point which we have also confirmed nu- 
merically. In other words, FFLO will extend all the way 
down to r s = in the phase diagram, which is contrary to 
what was found in [15] . Note that this is not the case when 
the interaction is screened, since this removes the singu- 
larity at k' = k, thus leaving an integrable singularity at 
Q = k' = 2k/r s . 



p-4 



Supersolidity in electron-hole bilayers with a large density imbalance 



Dilute gas of minority particles. To understand 
the implications of our single-minority-particle phase dia- 
gram for a dilute gas of minority particles (and, thus, the 
density-imbalanced electron- hole bilaycr), we must deter- 
mine the interaction between minority particles. Firstly, 
in the normal phase, the effective interaction V 22 be- 
tween two unbound minority particles within RPA satisfies 
V 22 = C/q + g q flig q c [25] and therefore is given by 
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(10) 



In the long- wavelength limit (q — > 0), this interaction is 
repulsive and dipolar (V 22 (r) ~ 1/r 3 ), and thus we expect 
a dilute gas of minority particles (with density n 2 1/d 2 
when d 3> 1) to form a weakly-interacting Fermi liquid. 
This is similar to what was found in Ref. |26] for the mi- 
nority spins in a strongly-polarised 2D electron gas. 

The interaction between well-separated excitons is also 
easily calculated by treating the excitons as static dipoles 
and summing up the inter- and intra-layer contributions: 
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The exciton-exciton interaction in the vacuum is clearly 
repulsive and dipolar at large distances, and we see from 
Eq. (fTTj) that this remains the case in the presence of a 
Fermi sea. This suggests that FFLO excitons in the regime 
of large density-imbalance will be thcrmodynamically sta- 
ble against phase separation, unlike in the cold-atom case. 
However, note that our approximation for V£ x is no longer 
valid when the exciton size becomes comparable to the dis- 
tance between minority particles, which is the case near 
the FFLO-N transition. There is also the possibility that 
the short-range exciton-exciton interactions are attractive, 
in which case there will be biexciton formation. How- 
ever, QMC calculations for the two-exciton problem show 
that biexcitons cannot form when d > 0.25 for mass ra- 
tio a = 4 [37] ■ Indeed, biexciton formation will be even 
more suppressed in our spin-polarised case due to Pauli 
exclusion. 

Now that we have an estimate for the exciton-exciton 
interaction, we can determine the structure of the FFLO 
phase at large density imbalance. For a low density of mi- 
nority particles, we can treat the excitons as simple bosons 
and define the mean- field complex order parameter ip(r), 
where |VK r )| 2 corresponds to the exciton density. Fluc- 
tuations will, of course, prevent true long-range order at 
finite temperature in 2D, but our excitonic phases can still 
possess quasi-long-range order. 

The effective low-energy thermodynamic potential F[ip] 
is similar to that used in weak crystallisation theory |28j : 



F\ 



r/ 2 i 



-^| 2 + 7 |(V 2 + 



I 2 ■ 

* mtii 



Here, fi is the chemical potential for the excitons, where 
fi < describes an empty 2nd layer, while the term 



7 > sets the exciton's dispersion minimum at momen- 
tum Qmin ■ The repulsive interactions are taken to be con- 
tact, with A = Vq x — the momentum dependence of V^ x 
should not affect our solution below provided Q m ind <C 1, 
i.e. the typical length scale of the spatial modulations is 
large. Indeed, this condition is always satisfied within 
the FFLO regions of Fig. [T] For the FFLO phase, we 
consider solutions of the form ip(r) — X)n a « e " 1 '* 1 ' w ith 
|q?i| = Qmin- Substituting this into F[ip] we get 
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where A 



and B = ^ n a„a_„. Minimising 



with respect to the amplitude a n we find that the low- 
est energy solution (F[ipo]/£l = — /i 2 /3A) requires that 
B* ai = \B\at 1} B*a 2 = -\B\a*_ 2 , M = M = v^A, 
while a„ = for n > 2. Such a state is thus composed of 
four q„, i.e. ±qi and ±q2, so that we have: 



^o(r) = e 1 



[cos(qi • r) - i cos(q 2 • r + 9 12 )} , (12) 



Note that the phases 9, Q\ 2 and the directions qi, q 2 are 
randomly chosen. Equation (|12j) corresponds to an exciton 
condensate with a 2D spatial modulation — a supersolid. 
Contrast this with the opposite limit of weak polarization 
in the BCS regime, where the favored state is believed to 
be a single cosine [3J. 

Discussion. — The supersolid is expected to enjoy a 
sizeable region of existence away from Wigner crystallisa- 
tion, but one possible issue in the large imbalance limit is 
the formation of three-body states, or trions. These are 
known to exist for all mass ratios in the limits d — > and 
kp — > [29]. However, we expect trions to disappear with 
increasing d, like in the biexciton case, since the intra-layer 
repulsion will eventually dominate over the inter-layer at- 
traction. Moreover, we find that the spin-polarised trion is 
barely bound at d — and so we expect trions to exist only 
for d <C 1 in Fig. Q] We speculate that a spin-polarised 
electron-hole bilaycr (as considered in this work) may be 
better for achieving the FFLO phase than an unpolarised 
one, since spin polarisation both suppresses trion forma- 
tion and enhances the inter-layer attraction, due to the 
reduction in screening in Eq. (|B]). 

Thus far, our calculations have been restricted to zero 
temperature and we have ignored the effects of thermal 
fluctuations. However, the exciton binding energy Eb pro- 
vides a temperature scale below which the exciton should 
be robust against thermal fluctuations. Assuming that the 
minority particles are electrons (a ~ 0.25) and using the 
parameters in GaAs (ao — 7 nm and Eq ~ 17 mcV), we 
find that Eb for FFLO excitons at d/ao ~ 0.5 is of the 
order of 5 K near the SF-FFLO transition. Thus, FFLO 
excitons should be experimentally accessible. 
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In order to access the FFLO phase itself, we require a 
sufficiently large exciton density, i.e. a sufficiently large 
minority-particle density n%, since we expect the criti- 
cal temperature to scale with n 2 . Our predictions should 
be valid at finite n 2 provided the exciton size is smaller 
than the spacing between excitons, i.e. we must have 
n-z/ni < r 2 Es/E . Thus, our strong-coupling theory of 
FFLO should apply for n 2 jn\ < 1 deep within the FFLO 
region, e.g. for d/ao ~ 0.5 and r s ~ 10, but it will break 
down near the FFLO-N transition for finite n<i. Here, we 
expect a more BCS-like version of FFLO, with the FFLO- 
N transition shifting to lower r s as n 2 is increased [15] . 

For an ordinary excitonic superfluid in 2D, there will be 
a Berezinskii-Kosterlitz-Thouless transition [301131] to the 
superfluid state, with transition temperature given bjQ 



Tbkt 



Eq n 2 



(1 + a) 2 ksr 2 n\ 



(13) 



It is unclear whether or not this is also true for the FFLO 
phase, since translational as well as gauge invariance has 
been spontaneously broken. However, on general grounds, 
we expect the FFLO transition temperature to have the 
same scaling as Tbkt, so we can use to it to obtain 
an estimate. For the typical parameters r s ~ 10 and 
n 2 /n\ ~ 0.2, we get an FFLO transition temperature of 
order lOOmK, which is smaller than the exciton binding 
energy (as expected), but still within reach experimentally 
The FFLO phase we predict can be observed experi- 
mentally via light scattering off of the spatial modulations. 
In addition, if the electrons and holes are allowed to re- 
combine, then a signature of the finite momentum pairing 
will appear in the angular emission of the photons |33j . 
where we expect peaks in the exciton photoluminescence 
to occur at large angles with respect to the plane normal, 
corresponding to ±qi and ±q 2 . 
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